Enhancement of electron spin coherence by optical preparation of nuclear spins 
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We study a large ensemble of nuclear spins interacting with a single electron spin in a quantum 
dot under optical excitation and photon detection. When a pair of applied laser fields satisfy 
two-photon resonance between the two ground electronic spin states, detection of light scattering 
from the intermediate exciton state acts as a weak quantum measurement of the effective magnetic 
(Overhauser) field due to the nuclear spins. If the spin were driven into a coherent population 
trapping state where no light scattering takes place, then the nuclear state would be projected into 
an eigenstate of the Overhauser field operator and electron decoherence due to nuclear spins would 
be suppressed: we show that this limit can be approached by adapting the laser frequencies when a 
photon is detected. We use a Lindblad equation to describe the time evolution of the driven system 
under photon emission and detection. Numerically, we find an increase of the electron coherence 
time from 5 ns to 500 ns after a preparation time of 10 microseconds. 



Introduction. Single electron spins localized in small 
artificial structures, such as semiconductor quantum dots 
(QDs), have become available and to a large extent con- 
trollable P, H H 0] • Of particular interest is the phase 
coherence of such electron spins as single quantum ob- 
jects, both from a fundamental physics point of view and 
because of their potential use as quantum bits (qubits) 
for quantum information processing [(J \J\ • 

A number of physical mechanisms that lead to the 
gradual reduction of the quantum phase coherence (de- 
coherence) of the electron spin have been analyzed @- 
It has been established experimentally and theoretically 
that for an electron in a GaAs QD, the predominant de- 
coherence mechanism is the hyperfine co upling t o the nu- 
clear spins in the host material H GjJ [U 111 111 El For 
an unpolarized ensemble of TV nuclei and an effective hy- 
perfine interaction energy A, the dephasing time in the 
presence of a weak magnetic field is T 2 * ~ 1/er ~ ^J~N /A 
where a is the width of the distribution of nuclear field 
values h z parallel to the field. In a typical GaAs QD 
with A ~ 90/xeV or A/ gf i B = 3.5T 0, the number 
of Ga and As nuclei (spin 7 = 3/2) is iV ~ 5 • 10 5 and 
T 2 * ~ 5 ns- this value is supported by the experimental 
evidence H 13 . The T* decay originates from nuclear 
ensemble averaging and can be prolonged by narrowing 
the nuclear spin distribution |l3j |. Another strategy is to 
polarize the nuclear spins Q , but this requires a polariza- 
tion close to 100% which is currently not available [l3j . 
Two schemes have been proposed to achieve a narrowing 
of the nuclear spin distribution, based on electron trans- 
port E3 an d gate-controlled electronic Rabi oscillations 

Here, we analyze an optical scheme for nuclear spin 
preparation that makes use of spin-flip two-photon (Ra- 
man) resonance in a driven three- level system (TLS), 
in analogy to electromagnetically induced transparency 
(EIT) in atomic systems . The lowest electronic states 
in a QD formed in a III-V semiconductor (e.g., GaAs) 
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FIG. 1: Three-level system. State 1 (2) is a spin-up (down) 
conduction-band (Ec) electron, with splitting g^BB to t + Sh z 
where Sh z is the z-component of the nuclear (Overhauser) 
field fluctuations. State 3 is a trion with J z i — 3/2. Two 
laser fields with frequencies lu p and ui c are applied near the 13 
and 23 resonances with detunings Ai,2. For <r+ circularly po- 
larized excitation (along z ), both transitions are allowed for 
0^=0 and transitions to the J z i — —3/2 states are forbidden. 
Inset: Structural axis z' , leading to a splitting in Ey and spin 
quantization axis z || B to t in Ec where cos# — z ■ z' < 1. 



that are optically active under cr+ circularly polarized 
excitation are a Zeeman-split single electron in the lo- 
calized conduction band (Ec) ground state and the neg- 
atively charged exciton (trion) \X), i.e., two electrons 
(spin up and down) plus one valence band heavy hole 
(hh) with angular momentum J z i = +3/2 (Fig. [TJ. The 
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J = 3/2 sector in the valence band is split into lh and hh 
states along the quantization axis z' determined by the 
direction of strong QD confinement. Here, we assume 
excitation from the hh (J z > — ±3/2) subband only. The 
quantization axis z in Eq is parallel to the total magnetic 
field Btot, and we assume that the axes z and z' enclose 
an angle 9 > 0. The spin up and down states in Ec are 
then |t) = |T), - cos(0)|T) z '+sin(0)U)^ and \[) = \[ 
) z = cos(0)||) z ' — sin(0)| T) z '- From now on, we drop 
the index z. Two circularly polarized (<J+) continuous- 
wave lasers at the frequencies oj p — ujx — — Ai and 
lu c = lux — u I — A 2 stimulate the transitions between ||) 
and \X) and between ||) and \X), while the trion with 
J z i = —3/2 is not excited. 

The mechanism leading to a narrowing of the nuclear 
field distribution v works as follows. The population of 
the excited state \X), and thus the probability of photon 
emission (scattering) is only non-zero away from the two- 
photon resonance 5 — where 5 = A x — A 2 . In the 
presence of the nuclear spins this resonance moves to 5 — 
Sh z where Sh z is the deviation of the nuclear field (along 
z) from its mean (h z ). The absence of photon emission 
during a waiting time t constitutes a weak measurement 
of the quantum operator Sh z . In the limit t — > oo, it 
becomes a strong measurement, projecting the nuclear 
state onto \Sh z — 0) (width a = 0), thus eliminating 
electron decoherence due to the fluctuating field Sh z . 

Model. The Hamiltonian for the TLS coupled to nuclei, 



H = H + Hu 



(1) 



HloxPx, with Sj — 



and P x = \X)(X\ = 



contains the three energy levels in Hq - 

0~i 

x o 

(001) T (001). The spin splitting is given as hui z — 
gfisBtot — IffMsB + (h)|, the sum of the external mag- 
netic and the mean nuclear fields. The nuclear (Over- 
hauser) field operator is h = S i=1 ^Ii) where A{ — 
aiWo| v E'(ri)| 2 , and ^(ji) denotes the electron wave func- 
tion at the position of the ith atomic nucleus and vq 
is the volume of the unit cell. The classical laser fields 
in the rotating wave approximation (RWA) are described 
by El H int = fi p e^*|X)(t| + n^'lAOm + h.c. The 
coupling of the electron spin to the quantum fluctua- 
tions of the Overhauser field h, will be described by the 
term Hm = — ^<5h ■ S, where <5h = h — (h). In the ro- 
tating frame = U(t)i&(t) with U(t) = e~ l ^*P T + 
e — c * n + where p T = | T)(T | and p L = we 

find H(t) = U{t) [H(t) + huipP^ + Tiuj c Pi\ U(tf, and, up 
to constant terms proportional to the unity matrix (we 
drop the tilde and use H for the Hamiltonian henceforth) , 



H(t) 




n,. 



-A 



— 6h z Yi z 
2 



(2) 



where A = Ai + A 2 . The hyperfine flip-flop terms Hj 



ing rapidly at the frequency cu p 



" c)N l /4 are oscillat- 
gp B B to t/n - 8 



and can be neglected in the RWA [19( , leading to a block- 
diagonal Hamiltonian H = diag(iJ 1; H 2 , ■ ■ . , Hjc), with 



Hi, 



5h h z +S 








-8h k z 







(3) 



where 8h z are the eigenvalues of the operator Sh z . The 
state of the TLS combined with the nuclear spins is de- 
scribed by the density matrix p, which we divide up into 
3-by-3 blocks pkk' ■ The density matrix evolves according 
to the generalized master equation ^3 



p = Cp=-[H,p]+Wp, 
in 



(4) 



with the Hamiltonian Eq. J5J and the dissipative term 

Wp = J2a = T,l T Xa( 2(T aXPO~X a ~ 0~XXP ~ pa X x)/2 + 

Y J (3=i,xlp( 2(T Pf3P (T Pl3 ~ a Pl3P ~ P cr /3/j)/2, where cry = 
(Xij ® 1 = acts on the TLS only. The rate Vxa 

describes the radiative decay of the exciton \X) into one 
of the single-electron states a = | J.), while 7^ is the 
pure dephasing rate of state /3 = \ \),\X) with respect to 
|t). Since H has block form, the master equations for the 
various blocks are decoupled, and have the closed form 

pkk' = (H k pkk' - pkk'H k f) + Wpkk'- (5) 
in 

The diagonal blocks obey the familiar Lindblad equation, 
pkk = CkPkk, C k = -i[H k ,p] + Wp. (6) 

Stationary state. We start with the factorized initial 
state po = Xo ® v o with arbitrary initial density matrices 
Xo and vq of the TLS and the nuclear ensemble, where 
Vq = J2kk' Vkk'\Sh z )(5h z I, and \Sh z ) are eigenstates of 
the Overhauser operator, Sh z \Sh k ) = 8h k z \8h k z ). As the 
off-diagonal elements of vq turn out to be irrelevant for 
the stationary state of the TLS, our analysis is valid both 
for a pure initial state (e.g., with Vkw = s/vkkVk'k 1 ) and 
for a mixed initial state (e.g., the completely mixed state, 
Vkk' oc Skk')- We assume a Gaussian distribution Vkk = 
(27r)" 1 / 2 cr^ 1 exp[-(<5/^) 2 /2cr 2 ], with an initial width a = 
(T = A/y/N, plotted as a solid line in Fig. Ufa). For 
our numerical calculations we choose A = 90 peV and 
N w 5 • 10 5 corresponding to cr ~ 0.13/ieV ~ 0.2 KT 
with r = 1 ns. We choose a sample of n -C (21 + 1) N 
states from the total Hilbert space (n ~ 4000). 

Due to the hyperfine coupling between the TLS and 
the nuclei, the two systems become entangled as the sta- 
tionary state is reached. The stationary density ma- 
trix p is independent of ^0 and has the general form 
P = J2kk' P^k' ® \Shg){Shg |. We derived an analytical 
expression for the 3-by-3 diagonal blocks pkk of p as a 
function of all parameters, including Shk- We find numer- 
ically that the off-diagonal blocks pkk' (k ^ k'), i.e., the 
nuclear coherences, decay exponentially and thus vanish 
in the stationary density matrix p. Therefore, p is block- 
diagonal, pkk' = S k k' Pkk- 
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FIG. 2: Conditional evolution of the nuclear spin distribution 
u(8h z ) — Vkk- (a) During the first period t\ of evolution with- 
out photon emission, the initial Gaussian distribution (solid 
line) develops a peak at the two-photon resonance (dashed 
line), (b) Change of v{5h z ) after emission at t\ (solid line), 
until before emission time t2 of the second photon (dashed 
line). The two-photon resonance 8 has been shifted to the 
position of the left maximum (adaptive technique). The de- 
pleted region around 5h z = is developed at ti. (c) Anal- 
ogous situation between tu and £12- (d) v(8h z ) is obtained 
after a total time of Wps. Inset: Magnification of peak in (d). 
The width of v(8h z ) is reduced by a factor of w 100 compared 
to the initial width in (a). The parameters used for this calcu- 
lation are S7 C = Q p = 0.2 ns -1 , A = 0, Fx] = Txj. = Ins - , 
and 74. = jx = 0.001 ns~ . 



Evolution of the observed system. In order to enhance 
the electron spin coherence, we aim at narrowing the 
nuclear spin distribution. In the case of a Gaussian dis- 
tribution, this amounts to decreasing the width a with 
respect to its initial value o~q, thus increasing the electron 
coherence time to ~ l/2o~. Ideally, we would perform a 
projective measurement P on the nuclear spins such that 
PpkkP oc S(8h z — S). This also bounds the off-diagonal 
elements of v = Tr\p, where Tta is the partial trace over 
the TLS, because \vkk' | < \JvkkVk'k' due to the positivity 
of v. The measurement P can be successively approxi- 
mated by monitoring the photon emission from the QD. 
The longer the period t during which no photon is emit- 
ted, the higher is the probability for 8h z to be at the 
two-photon resonance, 5h z = S. 

To describe the state of the system conditional on a 
certain measurement record, we use the conditional den- 
sity matrix p c . In the absence of photon emission, p c 
obeys Eq. © with Ck replaced by Ck — S, with the col- 
lapse operator S |2fj |. 



Pkk = (A - S ) Pkk, S P= Y] FxaO-aXPO-Xa- (7) 
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FIG. 3: Electron coherence function |(S + (£))|/|(S+(0))| vs. 
electronic precession time t calculated from u(5h z ) in Fig. [5] 
after emission of the nth photon (n = 1,6,..., 26). The initial 
decay is approximately Gaussian. 



We have numerically calculated p c conditional on the ab- 
sence of emitted photons for a given duration t. We plot 
the updated probability distribution v^k from v = TrAp c 
as a dashed line in Fig. Eta). We find that the a posteriori 
probability is concentrated around the two-photon reso- 
nance. Of course, this process will eventually be stopped 
by the emission of a photon. 

Photon emission. The stationary emission rate is |20| 



r cm = TrSp(t) = ry^(pfcfc) 



kk)XX^kk, 



(8) 



where T — Tx| + Fxi- The average number of photons 
emitted during time t is (N p h) = tT em and the a priori 
probability for iVph = is, according to Poissonian statis- 
tics, Pdark(t) = exp(— r em t). The waiting time distribu- 
tion for photon emissions is p W ait(i) = T"^ exp(— r cm t) 
with a mean waiting time (t) = The progressive 

narrowing of the Overhauser field distribution, Eqs. Q 
and JSJ, leads to a decreasing photon emission rate r em , 
and therefore to an increasing average waiting time t. We 
first assume that every emitted photon is detected and 
later generalize to the case of imperfect detection. 

Using Eq. JSJ, we find for the update rule of the 
nuclear density matrix upon photon emission v = 
Ti\Sp c /TiSp c . The Overhauser field distribution after 
the emission is 



Vkk{Pkk) 



kkjXX 



'kk 



(9) 



where v^u and (pkk)xx = (X\pkk\X) are taken before the 
emission. According to Eq. (JHJ), the population in the 
Overhauser field Sh z corresponding to the two photon 
resonance 5h z = S is depleted by the photon emission 
(Fig. 12b, solid line). However, the combined effect of 
the observed evolution without emission and the collapse 
leads to narrowing of the nuclear distribution. 
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FIG. 4: Characteristic time to of the initial Gaussian decay 
in |(S+(t))|/|(S+(0))| in Fig. [31 as a function of the optical 
preparation time t, averaged over 50 numerical runs (with 
error bars indicating the standard deviation). 

Adaptive technique. The stationary, isolated TLS at 
the two-photon resonance is in a dark state. However, 
the coupling to the nuclei introduces a nonzero probabil- 
ity for occupation of the \X) state and for emission of a 
photon. The detection of a photon provides information 
about 6h z and an adjustment of the EIT setup. Thus, 
photon emission does not necessarily signify a failed at- 
tempt to narrow the nuclear field distribution, but can 
be used as an input for the next weak measurement with 
adjusted frequencies of the driving lasers, oj' p = uj p + e/2 
and uj' c = uj c — e/2, so that the new two-photon resonance 
condition is Sh z = 8' where S' = 5 + e while A' = A. 
We choose e such that the new resonance with the Over- 
hauser field lies in one of the two maxima <5ti™ ax formed 
after the photon emission, see Fig. |2t>- This situation 
is formally described by Eq. with the substitution 
S — > 5 + <Sti™ ax . We note that the adaptive technique 
also works by changing only one of the laser frequencies. 
Right after the photon emission, the TLS is in one of the 
single electron states, ||) or Within a time 1/r, much 
faster than any nuclear time scale, the driven electronic 
system will reach the new stationary state. Then, the 
photon emission from the QD can again be monitored, 
leading to an enhanced nuclear population at the new res- 
onance (Fig. [2b, dashed line), thus further narrowing the 
nuclear distribution and further enhancing the electron 
spin coherence. This procedure can be repeated many 
times, leading to a nuclear distribution that is limited 
only by the width of the EIT resonance (Fig. |2t,d). 

Electron spin decoherence. The electron spin coher- 
ence is quantified using the expectation value of the rais- 
ing operator S+(t) in a state \x+) that is prepared per- 



pendicular to the total field B to t and is freely process- 
ing about the fluctuating nuclear field Sh z , (S+(t)) = 
(x+\S + (t)\x + ) = J2 k u kk {Sh k z \{x + \S+(t)\x + )\5h k z ). Us- 
ing (6h*>\S+(t)\6h*) = e ltSh *°*/ 2 S+(0)e- u5h >^ 2 , we ob- 
tain (S+(t)) — (h/2)J2k l 'kkexp(itdh z ), which we plot 
in Fig. [3 at various stages in an adaptive optical mea- 
surement scheme. We make a Gaussian fit, (S+(t)) oc 
exp(— t 2 /ip), for short times t and plot the coherence time 
to as a function of total waiting time in Fig. 0] This is 
the main result of our theoretical analysis: The repeated 
observation of the QD photon emission and adaptation 
of the laser frequencies lu c and u> p after each photon emis- 
sion leads to a pronounced enhancement of the electron 
coherence time, for the realistic parameters chosen, from 
to = 5 ns to rs 500 ns within a total observation time of 

10 (AS. 

Imperfect detectors. We cannot expect to have perfect 
photon detectors at our disposal, therefore we discuss 
here the case of a detector with efficiency e < 1. For an 
imperfect detector, Eq. Q becomes p kk = (C k — eS) p kk , 
reflecting that a photon is only detected with probability 
e. We have numerically analyzed the case of e = 10% 
(other parameters as above) and find to sa 460 ns after a 
preparation time of t = 50 /is. The requirement for longer 
preparation times is a consequence of the reduced photon 
detection rate eT em . We note that this is still orders of 
magnitude shorter than the time after which the nuclear 
spin diffuses, estimated to be around 0.01 s due to higher- 
order hyperfine flip-flop terms [Tt) . but possibly longer 
due to Knight shift gradient effects. Long time scales for 
the decay of nuclear polarization, on the order of seconds 
in the case of nuclear spins in contact with donors in 
GaAs, have also been seen experimentally |2l| . 

Conclusions. We find that it is possible to efficiently 
enhance the quantum phase coherence of an electron spin 
in a QD surrounded by a large ensemble of nuclear spins 
by a continuous weak measurement of the Overhauser 
field using optical excitation at a two-photon resonance 
of the TLS formed by the spin-split conduction band elec- 
tron and one of the trion states. An intriguing question is 
whether the electron spin coherence can be enhanced by 
a quantum Zeno type effect to the point where it is ulti- 
mately determined by spin-orbit interaction: since reser- 
voir correlation time of dominant electron spin decoher- 
ence due to flip-flop terms of the hyperfine interaction 
is ~ 1/is, this would most likely require high efficiency 
detection of the scattered photons. 
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